PROCESS AND APPARATUS FOR CORRECTION OF A 



RESAMPLER 

BACKGROUND OF THE INVENTION 

[0001] This application claims a priority from German Patent 
Application No. 100 36 703.8, filed July 27, 2000, and the contents 
of the application are incorporated herein by reference. 

[0001 A] This invention concerns a process and an apparatus for 
the correction of resamplers. 

[0002] Resamplers, that is sampling-rate converters are used to 
convert a sampled digital input signal having one input sampling 
rate into a sampled digital output signal with a different output 
sampling rate. With a process of this invention there is an input 
signal having a sampling rate, or frequency, which is larger by an 
arbitrary factor (not necessarily by a whole number, or integer) than 
a symbol frequency or a chip frequency. With a WCDMA-signal, 
every data symbol in a chip sequence is coded so that each symbol 
is of a plurality of chips. A binary change between two amplitude 
values can take place between the chips. The invention, however, is 
also suitable for other digital signals with then the term "chip 
frequency" being replaced by the term "symbol frequency" or 
"symbol rate". Upon translating the input sampling rate into the 
symbol or chip frequency, the problem arises that the relationship 
between the input sampling rate and the symbol or chip frequency 
is only approximately known, since a timing generator of the 
resampler is not identical with a timing generator of the input 
sampling rate, and thus a drift between the two oscillators is 



possible. Further, an absolute phase position of the input sampling 
rate is unknown. 

[0002 A] It is therefore an object of this invention to provide a 
process and an apparatus for the correction of a resampler with 
which a sampled input signal, which is subjected to an input 
sampling rate and which has a symbol or chip frequency that differs 
from that of the input sampling rate, is translated into a sampled 
output signal in which a sampling rate corresponds to the symbol or 
chip frequency by changing the input sampling rate by a resampling 
factor, with which a drift between the input sampling rate and the 
symbol or the chip frequency as well as an absolute, constant 
shifting of the phase of the input sampling rate are compensated. 

SUMMARY OF THE INVENTION 

[0003] According to principles of the process of this invention, a 
non-linear operation is performed on the input signal Sd, so that a 
spectral line arises at the symbol or chip frequency fc; the input 
signal Sd is spectrally shifted by the symbol or chip frequency fc; the 
phase cp of the shifted spectral line is recorded at the symbol or chip 
frequency fc as a function of the sample time points; and a 
resampling factor (resamp_fac) or a time shift of the output signal Sc 
is corrected by a time correction value (resamp_offset) based on a 
regression of the phase cp of the shifted spectral line at the symbol or 
chip frequency fc as a function of the sample time points. 

[0003A] An apparatus according to principles of this invention 
comprises: a non-linear operating element that subjects the input 
signal Sd to a non-linear operation so that a spectral line results at 
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the symbol or chip frequency f c ; a frequency shifter, which 
spectrally shifts the input signal Sd by the symbol or chip frequency 
fc; a phase recording, or determining, device that determines the 
phase cp of the shifted spectral line at one of the symbol and chip 
frequency fc as a function of the sampling time points; and a 
regression and correcting device that, on the basis of a regression of 
the phase cp of the shifted spectral line at the symbol or chip 
frequency f c corrects the resampling factor (resamp_fac) as a 
function of the sample time points and/ or time- wise shifts the 
output signal Sc by a time correction value (resamp_offset). 

[0004] According to the invention, the input signal is subjected 
to the non-linear operation, for example a squaring. The non-linear 
operation produces spectral lines of the natural frequencies of the 
input signal. In this regard, a spectral line is produced at the 
symbol or chip frequency. Because the input signal can change its 
status only at the end of the symbol, or for a symbol of a combined 
plurality of chips of a CDMA signal only at the end of a chip, the 
input signal is modulated by the symbol frequency or the chip 
frequency and the symbol or chip frequency can be created as the 
spectral line by the non-linear operation. A further recognition 
leading to the invention is that by shifting the spectrum of the input 
signal such that the symbol or chip frequency falls near (in an error- 
free ideal case exactly) the rate zero, a particularly uncomplicated 
evaluation of the spectral line results. By determining the phase of 
the thusly shifted spectral line, an absolute time offset, that was in 
the input signal, as well as a relative time offset which adds from 
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sampling interval to sampling interval, can be directly estimated by 
linear regression. 

[0004A] Before the phase is determined, or recorded, preferably a 
decimating of the sampling value is performed by intermediate 
sampling with prior band limits. A filter that is used for this 
preferably has a frequency response with zero positions at the 
simple symbol or chip frequency and at double the symbol or chip 
frequency. Because of the prior spectral shifting by the symbol or 
chip frequency, the firstly cited zero position falls on the DC portion 
and the doubled symbol or chip frequency falls on the spectral line 
of the mirrored symbol or chip frequency. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0005] Further benefits, characteristics and details of the 
invention are explained in more detail below using an embodiment 
shown in the drawings. The described and drawn features can be 
used individually or in preferred combinations in other 
embodiments of the invention. The foregoing and other objects, 
features and advantages of the invention will be apparent from the 
following more particular description of the preferred embodiment of 
the invention, as illustrated in the accompanying drawings in which 
reference characters refer to the same parts throughout the different 
views. The drawings are not necessarily to scale, emphasis instead 
being placed upon illustrating principles of the invention in a clear 
manner. 

[0006] Fig. 1 is a block diagram showing how a correction 

apparatus of this invention is used; 
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[0007] Fig. 2 is a detailed block diagram of the correction 
apparatus of this invention. 

[0008] Fig. 3 is a time diagram explaining a drift and a 
necessary correction of an input signal; 

[0009] Fig. 4 is a graph of a spectrum of the input signal; 

[0010] Fig. 5 is a graph of a squared spectrum of the input 
signal; 

[0011] Fig. 6 is a graph of a squared and shifted spectrum of 
the input signal; 

[0012] Fig. 7 is a constellation diagram of the input signal; 

[0013] Fig. 8 is a graph of a spectrum of the squared and 
shifted input signal after a decimation of the sampling values; 

[0014] Fig. 9 is a graph of the phase progression, or curve, of 
the signal represented in Fig. 8 as a function of the sampling values; 
and 

[0015] Fig. 10 is a graph of the frequency response of a filter 
used to decimate the sampling values. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0016] Fig. 1 shows a portion of a receiving apparatus 1 with 
which a correction apparatus 2 of this invention is used. 

[0016A] An analog input signal S A is converted into a digital, 
sampled input signal S D in an analog/ digital converter 3. Following, 
in a decimating device 4, in this embodiment, the sampling rate is 
decimated, or reduced, by a factor of 2 and, in a thereon following 
receiving and correcting filter 5, whose function in this invention is 
not of further interest here, it is filtered and thusly fed to a 
resampler (sampling rate converter) 6. 

[0017] Since the sampled input signal only in the middle of the 
chips, from which the symbols of the WCDMA signal are composed, 
is of interest, the sampling rate in the resampler 6 is reduced to the 
chip frequency f c . Insofar as a sampling rate f A of the input signal 
S D relative to the chip frequency f c is subject to no drift, the 
relationship between the chip frequency f c and the sampling rate f A 
of the input signal S D is 25.6/3.84=6.66 in the shown embodiment. 
Because there is a drift in the sampling rate f A relative to the chip 
frequency fc, however, an error estimation of the chip rate fc and the 
absolute time of phase shift, which is represented in Fig. 1 as "chip- 
timing" or "timing offset to the chip time points", is necessary. The 
correcting apparatus 2 of this invention, which first estimates this 
error and corrects it, serves this purpose. 

[0017A] In the embodiment shown in Fig. 1, the deviation from 
the chip-rate or chip frequency f c is transmitted to a multiplier 7. 
This transmitted correction factor is 1 if there is no error so that a 

6 



resampling factor, resampjfac, transmitted to the resampler 6 is the 
ideal relationship between the input sampling rate £a and the chip 
rate fc in this case. If an error appears, the correction factor 
provided to the multiplier 7 correspondingly deviates from 1. In 
addition, a time shift, or timing-offset, resamp_offset, is transmitted 
to the resampler 6 so as to be able to compensate for a constant 
time offset. 

[0018] In Fig. 3 a static time shift, or timing-offset s, and the 
time offset drift, or timing-drift As, is shown. The timing offset 8 in 
the example shown in Fig. 3 totals 0.5, that is a half of a chip 
interval Tc. This constant, static timing offset s is in this example 
overlaid by a timing-drift As = 0.1. This means that the total time 
offset, which is composed of the static timing offset e and the timing- 
drift As increases from chip interval to chip interval by 0.1, that is 
by 10% of a chip interval Tc. 

[0020] The correcting apparatus 2 of this invention 
compensates for the static timing offset e as well as also the timing 
drift As. The estimation and correcting apparatus 2 of this invention 
is shown in more detail in Fig. 2. The digital sampled input signal 
Sd is fed to a non-linear operating element 8 that subjects the input 
signal Sd to a non-linear operation. This non-linear operation can, 
for example, be the forming of summed squares in which an in- 
phase component I and a quadrature phase component Q of the 
base band signal Sd are respectively squared and then summed 
(P+Q2). In a following multiplier 9, the factor & i-k-2jr ■ fAI fC i s 
applied to the output signal of the non-linear operating element 8, 

which means that a frequency spectrum of the output signal of the 
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non-linear operating element 8 is spectrally shifted by the chip 
frequency fc- The multiplier 9 therefore operates as a frequency 
shifter. In a following decimating device 10, in the shown 
embodiment, decimation of the sampling value takes place in a 
relationship 1/256, that is a sub-sampling in a relationship of 
1/256. In order to avoid an aliasing, an appropriate band limiting 
was previously performed. 

[0021] In a thereon-following phase recording, or determining, 
device 1 1 the phase of the output signal of the decimating device 10 
is recorded as a function of the sampling value, or function of the 
time. In this regard, it is important that no jumps take place at area 
boundaries of, for example, +180° to -180°, rather that the phases 
at the area boundaries be continuously updated. This can take 
place in a known manner, for example, by ignoring the transfer of 
arithmetic and logic unit, that is, unwrap, which is illustrated by 
element 12. 

[0022] Finally, a linear regression of the phase as a function of 
the sampling value is performed. A resulting compensation straight 
line resulting therefrom can, for example, be determined by the 
method of the smallest sum of the resilient squares (Summe der 
Federquadrate) . As is further shown below in more detail, the static 
timing offset or timing error s can be gotten from an axis portion of 
the compensation straight line. The timing drift Ae can be estimated 
from the slope of the compensation straight line. The linear 
regression is illustrated by the element 13. 
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[0023] The signal is described in more detail below using Figs. 
4, 5, 6, 8 and 9 respectively after each of the individual processing 
steps. 

[0024] Fig. 4 shows a spectrum of the input signal So sampled 
with the sampling rate £a as a function of a standardized frequency. 
Fig. 5 shows the spectrum at the output of the nonlinear operating 
element 8, with here the non-linear operating element 8 having 
carried out the summed squares. In this regard, three spectral lines 
are shown to result from the non-linear operation. A first spectral 
line 14, at a zero frequency, originates in the direct voltage portion 
(DC portion), which arises from the summed-squares operation. A 
second spectral line 15 is further evaluated with the process of this 
invention. Assuming the timing drift As is zero, this spectral line 
lies exactly at the chip frequency fc. In addition to this described 
spectral line 15 at -fc, a mirrored spectral line 16 arises at +fc. 

[0025] Fig. 6 shows a spectrum at the output of the frequency 
shifter 9. This spectrum agrees with that of Fig. 5, however, it is 
shifted by the rate fc so that the spectral line 15 is exactly zero, if no 
error appears and the timing drift As is zero. If the timing drift As 
differs from zero, the spectral line 15 deviates from the frequency 
zero. A possibility for determining this frequency deviation could be 
interpolation of the maximum of the spectral line 15 to directly 
determine the middle frequency of the spectral line 15. This 
solution, however, has proven to be relatively expensive. Thus, 
instead, according to this invention, it is suggested that a time- 
domain evaluation be performed in which the phase be subjected to 
a linear regression as a function of the time, or as a function of the 
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sampling time points. Before this, however, the number of the 
sampling values (samples) are reduced, or decimated. The spectrum 
of the decimated sampled values is shown in Fig. 8. By employing a 
band limitation during decimation a noise/ audible spectrum is 
narrowed so that a signal/ noise relationship is greatly improved, as 
is clearly shown by a comparison between the noise-amplitudes in 
Fig. 6 and Fig. 8, relative to the amplitudes of the spectral lines 15. 
This brings about also a reduced noise amplitude in the phase 
characteristic curve in the time domain. 



[0026] In Fig. 9 the phase cp of the signal is represented in the 
time range as a function of sampling time points (samples) at the 
output of the decimating device 10, that is, in the area of the 
spectral line 15. The linear rise of the phase characteristic curve, 
which is overlaid by a static noise, versus time, can be recognized 
therein. The compensation straight line 17 can be created, for 
example, by a minimizing of the sums of the spaced squares, or 
another regression process. When this is done the axis cpo is a 
measurement for the static timing offset 8. The conversion 
calculation can be performed with the formula 

(i) 

The slope of the compensation straight lines 1 7 is a measurement of 
the timing drift Ae, wherein the estimated timing drift Ae according 
to this procedure can be calculated according to the formula 

A<p — 

As= f c 'dec_fac {2) 
In 



10 



In this regard, Acp is the slope of the compensation straight line 17 
per the example, £a is the sampling rate of the input signal Sd, fc is 
the chip rate, or the chip frequency, and dec_fac is the decimating 
factor with which the decimating device 10 decimates the sampling 
frequency fA (in the example the decjac = 256). For illustration 
purposes, the change in the phase is shown in Fig. 9 as being above 
50 sampling values, that is 50 • Acp. It should be kept in mind that, 
in this regard, the spacing of the sampling values (samples) in Fig. 9 
must be converted into the scale of the period Tc of the chip time 
points, because the decimation of the sampling rate comes to 
dec_fac/fA, and as is provided by formula (2). 

[0027] For controlling the resampler 6, the resampling factor, 
resamp_fac, and the time shift (timing offset), resamp_offset, as is 
shown in Fig. 1, are used. The conversions into these control values 
result from the formulas: 

resamp_fac= f£ r (3) and 

/*•(! + A*) 

resamp_offset= e-— (4) 

fc 

By the estimation process of this invention for the static timing 
offset s and the timing drift Ae, the control values can therefore be 
created which correct the sampling rate and the absolute phase 
position of the resampler 6. 

[0028] For further illustration, a status diagram of the input 
signal Sd is shown in Fig. 7 without the correction of this invention. 
In this regard, a frequency offset additionally arises so that the 
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phase indicator in the diagram rotates and, therefore, a time 
characteristic curve of the control can be more clearly seen. It can 
be recognized that in an area 18 a relatively high scattering appears 
that converges in an area 19 and diverges in an area 20. With the 
correction of this invention there is a concentration of this 
scattering on an ideal circular line that after correction of the 
frequency offset can be reduced in the status diagram to the ideal 
four status points. 

[0028A] For further illustration, Fig. 10 shows the frequency 
response of the filter used in the decimating device 10. It can be 
recognized, in this regard, that the frequency response at the chip 
frequency fc and the double chip frequency 2fc, that is at the 
spectral lines 14 and 16, has zero positions, so as to suppress an 
influence to the estimation result from the spectral lines 14 and 16. 
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